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STABIUTY AND CONTROL DERIVATIVES 


OF THE T-37B AIRPLANE 


Mary F. Shafer 
Flight Research Center 


INTRODUCTION 


Because of the continuing interest in flight simulation and handling qualities , 
there is a requirement for reliable estimates of the stability and control derivatives 
of most types of aircraft . In response to these requirements , the NASA Flight Re- 
search Center perfected a technique for minimizing the effort of determining the 
stability and control derivatives of aircraft from flight data (ref. 1) and developed 
a set of FORTRAN computer programs to implement the technique (ref. 2) . The 
method of derivative extraction is based on a modified maximum likelihood estimator 
that uses the Newton-Raphson algorithm to perform the required minimization . 

These computer programs are currently being used at the Flight Research Cen- 
ter to obtain stability and control derivatives for a wide variety of aircraft . Among 
the aircraft studied is the T-37B airplane, a small jet trainer. This report presents 
the estimates of the derivatives for the T-37B airplane determined by the modified 
maximum likelihood estimation technique from flight data. 

These flight data were selected from maneuvers performed in the course of a 
multiple purpose flight test program. As a result, the entire flight envelope was not 
studied in the flight test program. In some instances, the incremental effect of a con- 
figuration was studied instead of all possible configurations. 


SYMBOLS 


Cj rolling-moment coefficient 

C^ pitching-moment coefficient 


normal-force coefficient 



r 

a 


P 


6 

e 

6 


r 

Subscripts: 


normal-force coefficient for zero angle of attack and zero elevator 
deflection 

yawing-moment coefficient 

side -force coefficient 

roll rate, deg/ sec or rad/sec 
pitch rate, deg/sec or rad/sec 
yaw rate, deg/sec or rad/sec 

angle of attack with respect to body axes , deg or rad 
angle of sideslip , deg or rad 
aileron deflection , deg or rad 

elevator deflection , deg or rad 

rudder deflection , deg or rad 


p, q, r, a, p, 

6 ,’6 ,6 
a e r 


partial derivative with respect to the subscripted variable 


DESCRIPTION OF THE AIRPLANE AND INSTRUMENTATION 


The T-37B airplane (figs. 1 and 2) is a small two seat twin engine subsonic jet 
with a low wing and retractable landing gear.. The primary control surfaces are 
ailerons , elevators , and rudder . The airplane also has flaps and a speed brake , 
but the effect of the speed brake was not investigated in this study . Details and 
specifications for the airplane are given in reference 3 . 

Airspeed, altitude, and the pertinent stability and control quantities were among 
the data recorded . Angles of attack and sideslip were measured by vanes on a nose 
boom. Data were acquired by means of a pulse code modulation (PCM) system, 
which converts analog signals to digital format. Standard passive analog filters at 
40 hertz were applied to all the data signals . The digital data were recorded on 
magnetic tape and telemetered to a ground station for real time monitoring and re- 
cording. 


TEST PROCEDURE AND FLIGHT CONDITIONS 


Data were gathered from two flights in which a total of 166 stability and control 
maneuvers were performed. Of these maneuvers, 73 were performed primarily with 
elevator input, 51 were performed primarily with aileron input, and 42 were per- 
formed primarily with rudder input. The maneuvers were performed so that the 
linearity of the airplane model could be maintained . Thus . the airplane was flown at 
a stabilized flight condition before the pilot initiated the control input for the maneu- 
ver . The pilot inputs approximated doublets . The peak to peak amplitudes of the 
control inputs ranged from 5° to 15° for elevator deflection, from 20° to 30° for ai- 
leron deflection, and from 20° to 50° for rudder deflection. 

Four airplane configurations were investigated. The configurations were zero 
flaps, gear up; half flaps, gear up; full flaps, gear up; and zero flaps, gear down. 


METHOD OF ANALYSIS 


A maximum likelihood estimator method of analysis was used to determine a com- 
plete set of linear stability and control derivatives from the maneuvers performed in 
flight. This method is called the modified maximum likelihood estimator and is fully 
described in reference 1 . The method , sometimes called the Newton-Raphson method , 
is an iterative technique that minimizes the difference between the measured aircraft 
response and the computed aircraft response by adjusting the stability and control 
derivative values used in calculating the computed response. The Newton-Raphson 
algorithm was used to obtain the minimizations . The method can be modified to in- 
clude a priori information from previous calculations, flight tests, or wind-tunnel 
tests . This modification is made by including a penalty for adjusting the unknown 
stability and control derivatives away from the a priori values . If new information is 
contained in a flight maneuver, the estimate of the derivaiive is affected only slightly 
by the a priori information . If no new information is contained in a maneuver, how- 
ever , the a priori value results . A complete description of the computer program 
used for the derivative extraction and FORTRAN listings are given in reference 2 . 

In addition to giving estimates of the derivatives, this method of analysis pro- 
vides uncertainty levels for each derivative . The uncertainty levels are proportional 
to the approximation of the Cramer-Rao bounds described in reference 1 and are 
analogous to the standard deviations of the estimated derivatives. The larger the 
uncertainty level, the more uncertain the validity of the estimated value. The uncer- 
tainty levels obtained for a derivative from different maneuvers at the same flight 
condition can be compared to determine the most valid . Therefore , the uncertainty 
levels provide additional information about the validity of the estimate of the deriva- 
tive . 


RESULTS AND DISCUSSION 


Four maneuvers yielded completely unsatisfactory fits. The results of these ma- 
neuvers, two of which were performed with aileron inputs and two with r;idder in- 
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puts, are not presented. As a result, estimates from 162 maneuvers (98 percent of 
those flown) arc presented; 73 performed primarily with elevator input, 49 performed 
primarily with aileron input, and 40 performed primarily with rudder input. 

The zero flaps , gear up basic configuration is the most common in normal flight 
and is therefore presented as the basis for comparison in figures 3 to 6 . Many basic 
configuration maneuvers produced estimates with nearly identical values for angles 
of attack from 2° to 3°. Not all these nearly identical values are presented in fig- 
ures 4 and 6. 

The longitudinal stability and control derivative estimates are presented in fig- 
ures 3 and 4, and the lateral-directionf 1 derivative estimates are shown in figures 5 
and 6. The symbol shows the value of the derivative for each maneuver, and the 
vertical bar associated with each symbol represents the uncertainty level for that es- 
timate . The estimates of Cj. and 6 do not have uncertainty levels , since they 

®trim 

are calculated from other estimates . Not all control derivatives are plotted for all 
maneuvers, because control derivatives can be estimated only when that control var- 
ies during the maneuver , and for come maneuvers only one control varied . 


Longitudinal Derivatives 

The estimates of the longitudinal stability and control derivatives for the gear up 
and gear down configurations with zero flaps are shown in figure 3 . Gear position 
had virtually no effect on . The changes due to gear position in all the other 

derivative estimates are apparent in figure 3. 

The effects of the various flap settings with gear up are shown in figure 4. Flap 
position had virtually no effect on C . All other derivatives show more marked 

changes. The data indicate that there is little difference for most of the derivatives 
between the effects of half and full flaps . The zero, half, and full flap settings re- 
sult in significantly different estimates for the normal-force coefficient, How- 

ever, all three fairings have approximately the same slope when Cj^ is plotted 

against angle of attack . The coefficient C.. is calculated from C„ , C., , and 

N N Nc 

a 6 

c 

. The values of vary most with flap setting and cause most of the offsets in 



Lateral-Directional Derivatives 

The lateral-directional stability and control derivatives for the gear up and gear 
down configurations with zero flaps are presented in figure 5 . The derivatives C, , 

P 

C , C , C » C, , C„ , and are not markedly affected by gear posi- 

• I II II^ lo • O ^ O 
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tion . The estimates of the other derivatives are affected by gear position . 


The estimates for the various flap settings with gear up are presented in fig- 
ure 6. The derivatives C„ , C , C, , , and C„ are not significantly 

^ a r r a 

affected by the position of the flaps. Only C. shows much difference between 

^6 


a 

half and full flaps. For the other derivatives, as in the longitudinal cases, the 
amount of the deflection of the flaps does not have a notable effect . 


Fairings 

The fairings of most of the stability and control derivatives are based solely on 
the estimates and the uncertainty levels associated with each of the estimates. In 
figures 5(f) and 6(f) , the fairings of C. are not consistent with all of the data 

r 

shown . These fairings were determined from the quality of the comparison between 
the computed roll rate and the measured roll rate for the portion of the maneuver 
where the rudder was varying. The C. estimates near the fairings resulted from 

r 

maneuvers where the comparison between the responses was good. 


CONCLUDING REMARKS 


A complete set of the linear stability and control derivatives of the T-37B air- 
plane was determined with a modified maximum likelihood estimator. The deriva- 
tives were extracted from subsonic flight data from two flights for the longitudinal 
and lateral-directional modes. Four airplane configurations were investigated: zero 
flaps, gear up: half flaps, gear up; full flaps, gear up; and zero flaps, gear down 
Of the 166 maneuvers flown, 98 percent yielded satisfactory results. 

The data indicate that the amount of flap deflection has a significant effect on the 
magnitude of the stability and control derivatives , although half and full flaps 
caused similar changes in most of the derivatives . Some of the significant deriva- 
tives were affected by the position of the gear. 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, California 93523 
September 11, 1975 
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Figure 1. T-37B airplane 
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Figure 4. Longitudinal stability and control derivatives 
for gear up configurations. 
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Figure 5. Continued. 
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Fijfure 5. Continued. 
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II 


Rap setting 

o Zero Fairing for zero flaps 

□ Half Fairing for half and full flaps 

O Full 



Fi^re 6. Lateral- directional stability and control derivatives 
for gear up confi^rations. 
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